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Abstract —We consider the fairness in cooperative data ex¬ 
change (CDE) problem among a set of wireless clients. In this 
system, each client initially obtains a subset of the packets. They 
exchange packets in order to reconstruct the entire packet set. 
We study the problem of how to find a transmission strategy that 
distributes the communication load most evenly in all strategies 
that have the same sum-rate (the total number of transmissions) 
and achieve universal recovery (the situation when all clients 
recover the packet set). We formulate this problem by a discrete 
minimization problem and prove its M-convexity. We show that 
our results can also be proved by the submodularity of the 
feasible region shown in previous works and are closely related to 
the resource allocation problems under submodular constraints. 
To solve this problem, we propose to use a steepest descent 
algorithm (SDA) based on M-convexity. By varying the number 
of clients and packets, we compare SDA with a deterministic 
algorithm (DA) based on submodularity in terms of convergence 
performance and complexity. The results show that for the 
problem of finding the fairest and minimum sum-rate strategy 
for the CDE problem SDA is more efficient than DA when the 
number of clients is up to five. 

I. Introduction 

Due to the growing amount of data exchange over wireless 
networks and increasing number of mobile clients, the base- 
station-to-peer (B2P) links are severely overloaded. It is called 
the ‘last mile’ bottleneck problem in wireless transmissions. 
Cooperative peer-to-peer (P2P) communications is proposed 
for solving this problem. The idea is to allow mobile clients 
to exchange information with each other through P2P links in¬ 
stead of solely relying on the B2P transmissions. If the clients 
are geographically close to each other, the P2P transmissions 
could be more reliable than B2P ones. 

Consider the situation when a base station wants to deliver a 
set of packets to a group of clients. Due to the fading effects of 
wireless channels, after several broadcasts via B2P links, there 
may still exist some clients that do not obtain all the packets. 
However, the clients’ knowledge of the packet set may be 
complementary to each other. Therefore, instead of relying on 
retransmissions from the base station, the clients can broadcast 
combinations of the packets they know via P2P links so as 
to help the others recover the missing packets. For this kind 
of transmission method, there is a so-called cooperative data 
exchange (CDE) problem: how to find an efficient transmission 
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strategy that achieves the universal recovery , the situation 
when all clients recover the entire packet set. 

Finding the minimum-sum rate strategy, the transmission 
scheme for universal recovery with the minimized total num¬ 
ber of transmissions among clients, is the most commonly 
addressed CDE problem. It was introduced in m and studied 
in 0-0. There are also other optimization problems, e.g., 
finding the strategy that minimizes the weighted sum of 
transmission costs 0 - 13 . However, the solutions to most of 
these problems are usually not unique. On the other hand, in 
P2P communications, we wish to distribute the communication 
load as evenly as possible to prevent running out of one (or 
several) clients’ battery usage and also promote incentives for 
the clients to cooperate in a fair manner. Therefore, a natural 
question that follows is how to find the fairest transmission 
strategy in a solution set. 

In this paper, we study the problem of finding the fairest 
solution among the constant sum-rate strategy set, the set 
that contains all transmission schemes that achieve universal 
recovery and have the sum-rate equal to a constant number. 
We formulate this problem by a discrete minimization problem 
and prove its M-convexity: The feasible region is an M- 
convex set, and the objective function is M-convex. We 
show that the solution can be searched by a steepest descent 
algorithm (SDA) based on the optimality criterion of M- 
convex functions. We analyze the differences in convergence 
performance and complexity between SDA and a deterministic 
algorithm (DA) proposed in 0 based on the submodularity. 
By applying both SDA and DA to the problem of finding the 
fairest solution in the minimum sum-rate strategies, we show 
that SDA converges faster and involves less complexity than 
DA when the number of clients is up to five. 

A. Related Works 

The fairness in CDE problem has also been studied in 
0, 0. In 0, the best solution in terms of Jain’s fairness 
index is found by solving a convex minimization problem 
among those strategies that minimize the weighted sum of the 
transmission costs of clients. The CDE system in (7[ allows 
packet splitting, i.e., 0 considers an optimization problem 
over continuous variables (real numbers). On the contrary, 
fH studies the fairness when packet splitting is not allowed, 
where a discrete minimization problem is formulated, and the 
submodularity of this problem is proved and used to propose 
a DA algorithm. This paper studies the same problem as in 
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TABLE I 

EXAMPLES OF DEFINITIONS OF fj 



{Pi, P2, P3, P4, P5} {Pl,P2,P6} 



type of fairness 

rj/a 2 

Jain’s fairness fiTI 

- log r j 

proportional fairness 1121 

rj log rj 

uniform fairness (8) 


Fig. 1. An example of CDE system: There are three clients that want to 
obtain six packets. The has-sets are Hi = {pi, P2, P3, P4, P5}, H 2 = 
{pi, p 2 , P 6 } and H 3 = {P 3 , P 4 , P 6 }• 


mi. The difference is that we prove the discrete convexity 
(M-convexity) of it, based on which an SDA algorithm is 
presented. In Section [TV] we show that the M-convexity can 
also be proved by the results derived in H] by the one-to- 
one correspondence between M-convex sets and submodular 
functions. In Section [VI we show the performance of SDA by 
comparing it to DA by examples. 

II. System Model and Problem Statement 

Let V = {pi,...,pat} be the packet set containing N 
linearly independent packets. Each packet p^ belongs to a 
field F q . The system contains K geographically close clients. 
Define the client set as /C = {1,..., K}. Each client j G JC 
initially obtains Hj C V. Here, Hj is called the has-set of 
client j. We also denote H C j — V\ Hj as the packet set that 
is missing at client j. The clients are assumed to collectively 
know the the packet set, i.e., VJj^jcHj = V. The P2P links 
between clients are error-free, i.e., any information broadcast 
by client j can be heard losslessly by client j' for all j' such 
that j' G JC\{j}. The clients broadcast linear combinations 
of the packets in their has-sets in order to help each other 
to recover the entire packet set V. For example, in the CDE 
system in Fig. HE client 1 broadcasting pi + P 3 helps client 2 
recover P 3 and client 3 recover pi, and client 2 broadcasting 
pi + p 6 helps client 1 recover p 6 and client 3 recover pi. 

Let r = (7*1, ..., vk) be a transmission strategy, where rj 
denotes the total number of linear combinations transmitted 
by client j. We call J2je)C r j sum-rate of strategy r. Let 
a be the minimum sum-rate that allows universal recovery. 
Denote a as an integer constant and assume that a > a. 
Consider the CDE problem when the sum-rate is constrained 
to a budget a , i.e., the problem of finding a transmission 
strategy that achieves universal recovery and has a sum-rate 
equal to constant a. The solution to this problem is not unique 
in general. For example, the CDE system in Fig. [I] has a = 4. 
Let a ==fc a «= 4 and consider two transmission schemes: one 
is that client 1 broadcasts pi + P 3 , P 2 + P 4 and P 5 and client 
3 broadcasts p 6 ; The other is that client 1 broadcasts P 2 + P 4 
and P 5 , client 2 broadcasts pi +p 6 and client 3 broadcast P 3 . 
The transmission strategies associated with the two schemes 
are (3,0,1) and (2,1,1), respectively. Both strategies achieves 
universal recovery and have sum-rate equal to 4. However, 

^he CDE system in Fig. [I] is also the example considered in |8 |. We 
use this system because we will compare our work with the results in ( 8 ) in 
Section HyiandlYl 


strategy (3, 0,1) is not as good as (2,1,1) in terms of fairness: 
In strategy (3, 0,1), the energy consumption at client 1 is high 
while client 2 is a free-rider. So, there arises a problem of how 
to find a fairest solution in the constant sum-rate transmission 
strategy set. 


III. Discrete Convex Minimization 

A. Discrete Minimization Problem 

It is proved in (9[, iflOll that the necessary and sufficient 
condition for a transmission strategy r to achieve universal 
recovery is that X^, e <s r j > I Ojeic\s I f° r ^ such that 
S C JC. Denote 

r (' S ) = L! r i- (1) 

jes 

Let 7 Z a be the set that contains all strategies that allow 
universal recovery and have sum-rate equal to a. We can 
describe 7 Z a as 


TZa = {r e Nq : r(<S) > | f] U) 


je>c\s 


,V<5 C JC, 

LA =Q! }- (2) 


je >c 


Note, for all a > a, set 7 Z a is nonempty, and when a = a, 
the problem under consideration is to find the fairest solution 
in the minimum sum-rate strategy set 7 Z&. 

Let F a : 7 Z a Gibe the fairness measurement functions. 
We assume that F a is a discrete separable convex function in 
1Z a , i.e., F a is in the form of 


F a (r) 


Xje/e/iG) r G n <y 
+00 r ^ lZ a 


( 3 ) 


where f j is convex in r 7 for all j. In TABLE U we show some 
examples of fj based on different types of fairness indices. 
In this paper, we use the uniform fairness definition fj (rj ) = 
rj log rj c| The fairest strategy that achieves universal recovery 
and has sum-rate equal to a can be searched by solving the 
minimization problem 


minF a (r). (4) 

r 

Denote dom z / = {x G 7L K : — 00 < /(x) < + 00 } the 
effective domain of function /. dom z F a = 7 Z a , i.e., dom z F a 
denotes the feasible region of 0. 

2 Although we just consider fj(rj) = rj log rj, it should be clear that the 
results derived in the paper is applicable to other definitions in TABLE U 
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B. M-convexity 

In this section , we show the M-convexity of 0]). We first 
clarify some definitions as follows. 

Definition 3.1 (M-convex set K73\l): A set B C Z> K is M- 
convex if it satisfies the following exchange axiom: 

B-EXC_[Z]: For all x, y G B and u G supp + (x — y), there 
exists v G supp - (x — y) such that x — e u + e v g8. 

In B-EXC_[Z], supp+(x) and supp (x) are the positive and 
negative supports of x, respectively, and defined as 

supp+(x) = {j: Xj > 0, j € 1C}, 
supp _ (x) = {j: xj < 0,j € 1C}, 

where xj is the jth entry of x. ej is the unit vector with all 
entries being 0 except the jth entry being 1 . 

Definition 3.2 (M-convex function 4731/): A function 
/: Ij K ^ R+ is M-convex if dom z / 7 ^ 0 is an M-convex 
set and for all x, y G dom z f and u G supp + (x — y) there 
exists v G supp - (x — y) such that 

/(x) + /(y) > /(x - e u + e v ) + /(y + e w - e v ). (5) 


Remark 3.3: M-convex set is a class of discrete convex 
sets based on B-EXC_[Z]. Let B C R K be the convex hull 
of B. If B is an M-convex set, B = B D 7L K , i.e., all integer 
points contained in B constitutes B. Alternatively speaking, B 
describes a hole-free region in 7L K [031 . M-convex function 
is a class of discrete convex functions defined based on B- 
EXC.[Z]. We will show examples of M-convex sets and 
functions in Example 13.91 

In the following context, we prove the M-convexity of (0]) 
by showing the M-convexity of !Z a and F a . We start the proof 
by showing a property of the tight set as follows. 


njc 

ijeK\s n j 


We call S C JC a tight set of r if r(<S) = P| ? . 
tight set has the following property. 

Lemma 3.4: Let r G lZ a and y C 1C such that XC\ y 
0 and X U y 7 ^ JC. if X and y are tight sets, then X ny is a 
tight set. 

Proof: Since r(X) + r(y) = r(X D 30 + r(X U 3), we 

have 


c(X n y) = r(X) + r (3) - r(X u y) 


3 n 

jeic\x 

n 




n n ^ 


jeic\y 


jeic\xuy 


< 


H c „ 


jeic\xny 


( 6 ) 


But, r{x ny)> 

n 


n 


\jeK\xny Wj 


since r e 7 Z 0 


Therefore, 

H y) = C\jeK\xc\y Hj . Note, the last inequality in © 
is because of Corollary IA.2I in Appendix ■ 

Theorem 3.5: 7Z a is an M-convex set. 

Proof: We use an approach similar to the proof of 
Proposition 4.14 in m. Assume that the exchange axiom 
B-EXC_[Z] does not hold in lZ a , i.e., for some x, y G lZ a , 
there exists u G supp + (x — y) such that x — e u + e v £ lZ a for 
all v G supp - (x —y). Since the sum-rate of x — e u + e v equals 
a always, the only situation that could make x— e u +e v £ lZ a 


is that for each v G supp (x — y), there exist S v C JC such 
that ueS v ,v(£S v and 


X 6^ + 


x(5 w ) - 1 < 


n 

j£JC\S v 


(7) 


But, x(S v ) > 


rije/c\s v 


n je/c\s v Wj 


Therefore, x(5 v ) 


, i.e., S v is a tight set for all v G supp (x —y). 
Consider the set Z = n V Gsu PP -(x-y) S v- Since S v D S v > / 0 
and 7 ^ JC for all v, v' G supp - (x —y). By Lemma [3~4l 

Z is a tight set, i.e., x(Z) = C\jejc\z • Also, since u G Z 
and v £ Z for all v G supp - (x — y), 


y (z) < x(z) = 


n w, 

jeK\z 


( 8 ) 


® contradicts the condition that y G lZ a . Therefore, for all 
u G supp + (x — y), there must exist v G supp - (x — y) such 
that x — e u + e v G !Z a , i.e., B-EXC_[Z] is satisfied in lZ a . 
By Definition 13.11 lZ a is an M-convex set. ■ 

Theorem 3.6: F a (r) is an M-convex function. 

Proof: Since each discrete separable convex function is 
also M-convex ED, F a (r) is an M-convex function. Also, 
dom z F a = lZ a is nonempty and M-convex. By Definition 13.21 
F a ( r) is an M-convex function. ■ 

Corollary 3. 7: © is an M-convex minimization problem. 

Proof: This is a direct result of Theorem 13.61 ■ 

Remark 3.8: Similar to continuous convex minimization 
problems, for M-convex minimization problems, local opti¬ 
mality guarantees global optimality, which is based on the 
optimality criterion ED 


F a (r) < F a (r - e u + e„),Vr G dom z F a . (9) 

Example 3.9: Consider the CDE system in Fig. [T] when a = 
4. We have 

n 4 = {(2,1,1), (3,0,1), (3,1,0)}. (10) 

It can be shown that IZ 4 satisfies B-EXC_[Z]. For example, 
consider x = (2,1,1) and y = (3, 0,1). supp + (x — y) = {2} 
and supp - (x — y) = {1}. x — e 2 + ei = (3,0,1) G IZ 4 . It can 
be checked that this property applies to all x, y G IZ 4 . When 
a = 5, we have 


n b = {(1,2,2), (2,1,2), (2, 2,1), (3,0,2), 

(3,1,1), (3,2,0), (4,0,1), (4,1,0)}, (11) 

where B-EXC_[Z] also holds. In Fig. [2 we show the set 7 ^ 4 , 
1Z 5 and their convex hulls IZ 4 and IZ 5 . It can be seen that IZ 4 
and 1Z 5 lie on the planes r 3 ~ 4 and r j = 

respectively, and IZ 4 and 1 Z$ are hole-free, i.e., all the integer 
points lie on IZ 4 and 7Z$ belong to IZ 4 and 1 Z$, respectively. 
Consider the function F^r) in Example 13.91 We have 


^4(r) 


'1.3963 r =(2,1,1) 

<3.2958 r= (3,0,1) or r= (3,1,0) . 
+oo otherwise 


( 12 ) 


Therefore, arginin r F 4 (r) = {(2,1,1)}. Similarly, we can 
show that argmin r F 5 (r) = {(1,2, 2), (2,1, 2), (2, 2,1)}. It 
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• U 5 



Fig. 2. IZ 4 and 7Z$ of the CDE problem in Fig. |T] IZ 4 and TI 5 are the 
convex hul l of TI 4 and TZs , respectively. It can be seen that IZ 4 = 72 4 D Z K 
and 7^.5 = 7^-5 n Z K . 



Fig. 3. The submodular polyhedron P(g 4 ) and submodular base polyhedron 
B{g 4 ) for the CDE problem in Fig. |T| g 4 is given in fl 6 l . It can be seen that 
is exactly IZ 4 , the convex hull of IZ 4 and P 4 contains all the integer 
points in P(g 4 ), i.e., ZI 4 = B(< 7 4 ) H Z K . 


can be verified that F 4 and P 5 satisfy the optimality criterion 
in Remark [3751 

IV. Relationship with Existing Works 

In 13 , H), the properties of the set lZ a has been studied. 
They both show that 7 Z a is related to a submodular set 
function. In this section, we show that the Theorem 13.51 can be 
proved by the results in ei, m and ( 0 ]) is in fact a constrained 
resource allocation problem. 


According to Theorem 2.6 in fPflL B(g a ) = B(g a ), and g is 
submodular. Alternatively speaking, 7 Z a = B(g a ) D Z K , i.e., 
!Z a can be fully determined by submodular set function g a . 
In the following context, we discuss the relationship between 
the submodularity of g a and the M-convexity of 7 Z a . 

Theorem 4.3 (one-to-one correspondence mi): The class 
of submodular set functions / and the class of M-convex 
sets B are in one-to-one correspondence through the mutually 
inverse mappings: 


A. M-convex set and submodularity 

We first clarify the associated definitions as follows. 

Definition 4.1 (submodular set function KT4\l): Let 2 ^ be 
the power set (the set of all subsets) of JC, f: 2 K ha M + 
is submodular if for all X, y C K 

f(x) + f(y) > f{x ny) + f(x u y). (13) 

Definition 4.2 (polyhedron and base polyhedron HBl): 

For a function /: 2 k ^ M + , the polyhedron P(f) and base 
polyhedron B(f) are defined as 

P(f)~-{reM%:r(S)<f(S),SClC}, 


f ^ B: B= B{f)nZ K , 

B 4 /: f(S) = sup{r(<S): r e B}, VS C/C. 

In Theorem 14.31 the base polyhedron B(f) is exactly the 
convex hull of B B3. Therefore, the interpretation of mapping 
/ ha S is: For every submodular function, its base polyhedron 
determines the convex hull of an M-convex set. Based on 
Theorem l43l Theorem [331 in Section HIl-B I can also be proved 
by the results in @, ( 8 ). Since lZ a = B(g a )C\Z K , due to the 
submodularity of g a , !Z a is M-convex. 

Example 4.4: In the CDE system in Fig. |T] The function 
g a determined by (fl5l) when a = 4 is m 


If / is submodular, P(f) and B(f) are submodular polyhe¬ 
dron and submodualr base polyhedron, respectively. 

In ©, m, it was shown that lZ a = B(g a ) D 7L K , where 


9a{S) 


f° 


[a- 

fW\ S u c j 


<S = 0 
otherwise 


was a crossing 


submodular function 0 Let 


(14) 


g a (S) = min( {V» : i e 1 = {1, 2 ,... }} 

iez 


is a partition of 



(15) 


3 a crossing submodular function / satisfies the inequality GU for all 
X,y<ZK such that X n y ^ 0, X - y + 0, y - X ± 0 and X U y + K. 


54(0) = 0, <? 4 ({1}) = 3, r?4({2}) = 1,5 4 ({3}) = 1, 
$ 4 ({1,2}) = 4, p 4 ({l,3}) = 4, p 4 ({2,3}) = 2, 

5 4 ({ 1 ) 2 , 3 }) = 4 . ( 16 ) 

We show the polyhedron P(g 4 ) and base polyhedron P(# 4 ) 
in Fig. 0 It can be seen that B(gfi) is the intersection between 
P(g 4 ) and plane J2 3 eJC T? = 4- By comparing Fig.[3]to Fig. [2 
it can be seen that P(g 4 ) is in fact P 4 , the convex hull of the 
M-convex set P 4 , and P 4 = P(g 4 ) D 7L K . 

In fact, this paper and 0, ( 8 ) present two different ways of 
proving the M-convexity of lZ a : We prove the M-convexity 
of lZ a by definition; The authors in [ 0 , [81 express lZ a by 
a submodular base polyhedron, which implicitly proves the 
M-convexity of lZ a by Theorem 14.31 
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Algorithm 1: Steepest Descent Algorithm (SDA) fl3l 

k = 0 ; 

Find a r (0) G U a \ 

repeat 

Find the descent direction — e u * + e v * by determining 
(u*,v*) such that 

G argmin{F a (r (/c) - e u + e v ): 

+ e u - e v G 1Z a }; 

r (fc+i) _ r (fc> _ 6u , +e v *; 
until F a ( r (/e+1) ) > F a ( r (fc) ); 

return minimizer r* = r^; 


B. M-convexity of Resource Allocation Problems 

Since the concepts of M-convex set and submodular base 
polyhedron are exchangeable based on the one-to-one corre¬ 
spondence in Theorem 14.31 the results in Section IIII-BI also 
implicitly indicate the relationship between problem © and 
resource allocation problems under submodular constrtains. 
Based on the results in Section HY-A1 we can rewrite problem 
© as 

minji^r): r e £(&*)}. (17) 

This is called the sparable convex resource allocation problem 
under submodular constrains ca (since the feasible region 
can be described by a submodular base polyhedron). The M- 
convexity of this problem has been proved in ||T3l . lfT 6 l . In 
m-m , various algorithms based on M-convexity have been 
developed for solving problem <H3. 

V. Steepest Descent Algorithm 

One of the advantages of solving an M-convexity mini¬ 
mization problem is that efficient methods can be derived. 
This section considers one of the simplest methods, the 
steepest descent algorithm (SDA). We analyze the convergence 
performance and time complexity of SDA by comparing it to 
the deterministic algorithm (DA) proposed in ( 8 ). 

The SDA shown in Algorithm 1 is directly devised based 
on the optimality criterion ©. It starts with r^°\ an arbitrary 
point in 7 Z a , and moves along the descent direction in each 
iteration. The authors in 0 also studied the fairness in CDE 
problem, where algorithm DA as shown in Algorithm 2 was 
proposed for solving problem ©. SDA differs from DA in the 
following aspects: 

(a) Running SDA requires the value of a and a starting point 
r(°) in 1Z Q ; Running DA only requires the value of a. 

(b) For SDA, G 7^ a ; For DA, rW G P{g a )- Since = 
B(g a ), SDA searches the minimizer in the submodular 
base polyhedron, while DA searches the minimizer in the 
submodular polyhedron. 

(c) The number of iterations in SDA is bounded by —— ~ r ^ 
E3); The number of iterations in DA is exactly a. 

Aspect (a) implies that if only a is known DA can be 
run independently while SDA requires the running of another 
algorithm to find i*(°) G lZ a . However, finding r^ 0 ) is not 


Algorithm 2: Deterministic Algorithm (DA) 0 

r (0) = 0; 

for k = 1 to a do 

Find j* such that 

j* € arg mm{fj (r j +1) — fj (r j) : j G K,r+ej G P(gf)}] 
r (fc+i) r (fe) + ej *; 

endfor 


difficult. It can be accomplished by the randomized algorithm 
proposed in 0. If the problem is to find the fairest solution 
in the minimum sum-rate strategies, i.e., a = a, there are 
many existing methods that find the value of a, and most of 
them returns a with a strategy i*(°) G 7 Z&, e.g., the divide- 
and-conquer algorithm in (4), the randomized algorithms in 
0, 0. It should be also noted that if SDA and DA are 
applied after obtaining a and the initial point G H a by 
the methods in 0-0 the knowledge of will be discarded 
in DA. Aspect (c) is in fact a result of aspect (b). It implies 
that SDA converges faster than DA if 

Li(a) = max{||x-y||i: x,y G H a }, (18) 

the Zi-size of 7 Z a , is smaller than a. 

In SDA, checking whether —e u +e v G H a is equivalent 
to checking whether rW — e u + e v G P(g a )• Therefore, we 
can use the same feasibility checking algorithm as in 0. The 
method is to run a submodular function minimization (SFM) 
algorithmic 

Example 5.1: Consider problem © when a = 4 for the 
CDE system in Fig. [I] By applying SDA algorithm, we 
get the estimation sequence {r^} = {(3,1, 0), (2,1,1)}, 
where the starting point i*(°) = (3,1,0) is obtained 

by running the algorithm in 0; By applying DA al¬ 
gorithm, we get the estimation sequence {r^} = 

{(0,0,0), (1,0,0), (1,1,0), (1,1,1), (2,1,1)}. It can be seen 
that SDA converges faster than DA. We plot the searching 
paths indicated by sequence {r^} for both SDA and DA in 
Fig. 01 It clearly shows that SDA searches the minimizers in 
the submodular base polyhedron B(gf), while DA searches 
the minimizer in the submodular polyhedron P{gf). We then 
run SDA and DA for a = 5 and a = 6. We plot ||^ ( (0) ) ~^|| , the 
normalized errors, in Fig. 0 It shows that SDA requires less 
number of iterations than DA before reaching the minimizer 
r*. 

Although SDA converges faster than DA, the complexity 
of SDA could be higher than that of DA. Fet ( denote the 
complexity of running SFM algorithm for the feasibility check. 
We assume that ( includes the complexity of the measurement 
of the objective function if the results of the feasibility check is 
true. For SDA, the complexity is bounded by 0(K 2 • £ • L\ (a)) 
urn Note, 0(K 2 • ( • L\{a)) is the maximum complexity of 

4 In both 0, 0, it was shown that checking whether r G was 
equivalent to a submodular function minimization problem. This paper applies 
the same feasibility checking algorithm as in 0, 0. For more details, we 
refer the reader to Algorithms IV. 3 and IV. 4 in 0. 
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Fig. 4. The searching paths of SDA and DA for problem @ with a = 4 in 
the CDE system in Fig. |T] The number of iterations is 1 for SDA and 4 for DA. 
For SDA, rW e B(g 4 )nZ K , i.e., r<» e fc 4 ; For DA, rW G P(g 4 )nZ K . 



Fig. 5. Convergence performance of SDA and DA for problem in the 
CDE system in Fig. \l\ || » the normalized error, when a = 4, 5 and 

6. It can be seen that the number of iterations in DA equals to the value of 
a. The number of iterations in SDA depends on the distance between the 
starting point r(°) and minimizer r*. It is less than DA for all three values 
of a. 

SDA for two reasons: one is that L\{a) is the maximum l\- 
norm of 7Z a ; The other is that we only need to run SFM 
algorithm (to check the feasibility) for those (u,v) such that 
r W — e u + e v > 0. For DA, the complexity is exactly 0(K • 
( • a). Therefore, the main difference in complexity between 
SDA and DA is the number of runs of SFM algorithm. Since 
Li(a) grows with a , the computation load of SDA may be 
heavier than that of DA when L\ (a) is comparable to or larger 
than a. 

Consider problem © when a = a , i.e., the problem of 
finding the fairest solution in the minimum sum-rate strategies. 
In this case, L\(a) is minimum, and a is proportional to N, 
the number of packets d). Therefore, the computation loads 
of SDA and DA grow with both K and N, and the complexity 
growth in SDA is larger than DA. 

Example 5.2: We run an experiment to compare SDA and 
DA in terms of convergence performance and complexity 
when they are applied to search the fairest solution in the 
minimum sum-rate strategies. We vary K, the number of 
clients, from 3 to 10 and N, the number of packets, from 
6 to 30. For each combination of K and N, we repeat the 



(a) Convergence performance of SDA and DA: the average number of 
iterations before reaching the minimizer r* 



(b) Time complexity of SDA and DA: the average number of runs of 
SFM algorithm 

Fig. 6. Comparisons of SDA and DA when applied to problem a with 
a = a, the problem of finding the fairest solution in the minimum sum-rate 
strategies. K, the number of clients is varying from 3 to 10, and N, the 
number of packets, is varying from 6 to 30. 


following steps by 20 times. 

• Randomly generate the has-set ELj for each client subject 
to the condition U jelc'Hj — 'P- 

• Run the randomized algorithm in © to obtain the mini¬ 
mum sum-rate a and the starting point E 

• Apply SDA and DA to search the fairest transmission 
strategy r*@ 

The number of iterations and complexity averaged over repe¬ 
titions are recorded and shown in Fig. [6j In Fig.^a), it can be 
seen that SDA always converges faster than DA and there is a 
clear growth in the number of iterations in DA with increasing 
According to Fig. Ob), the complexity of SDA is lower 
than or close to that of DA when K is lower than 6. It can be 
seen that the complexity of SDA is much higher than that of 
DA when K and N are large. The results in Fig. 0 shows that 
SDA is more efficient than DA when the number of clients is 
up to 5. 

5 The randomized algorithm in 0 returns a with a strategy r(°) E 7^. In 
SDA, both a and r(°) are used. But, in DA, only a is used. 

6 It also implies that a grows with K more drastically than L± (a). 
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VI. Conclusion and Future Work 


We formulated a discrete minimization problem for finding 
the fairest solution in the constant sum-rate strategies that 
achieved universal recovery in CDE system. We proved the 
M-convexity of this problem and presented an SDA algorithm 
for searching the minimizer. We discussed the relationship 
between our work, the results in m and resource allocation 
problems with submodular constraints. By a comparison in 
convergence performance and complexity between SDA and 
DA, we showed that SDA was more efficient than DA when 
it applied to the problem of finding the fairest solution in the 
minimum sum-rate strategies among a small number of clients. 

As part of the conclusion, we briefly discuss how the results 
in this paper can guide the research work on CDE problems 
in the future. One can study the M-convexity of other opti¬ 
mization problems over the feasible region 7 Z a . For example, 
the weighted sum transmission cost minimization problem 
min{w T r: r G where w is a weight vector. Since w T r is 
separable convex, the M-convexity of this problem is directly 
proved by Theorems 13.51 and 13.61 in Section [HE UH For solving 
M-convex optimization problem, one can consider algorithms 
other than SDA, e.g., the algorithms proposed in lfT71 - llT9li FI 
The complexity in these algorithms could be lower than SDA, 
i.e., they could be more efficient than both SDA and DA. 


Appendix A 


Lemma A.l: For Hj C V and X, y c JC such that XfAy / 
0 and * U y ± /C, 


In 

jex 


Hi 


n«s 

jey 


< 


n 

jexuy 


+ 


n 

jexny 


Proof: Recall that f |jexF ^ fWn yF and 
Hiey C Ojexny because Xny C X and Xny C y, 
respectively. We have 


|(pl^i)u(pl^i)| 

jex jey 

<|( n ^) u ( n u i) 

jexny jexny 

=|( n 

jexny 


(19) 


1 F a ( r) and w T r are both separable convex in r. But, separable convex 
function is just a special case of M-convex function. For the problem where 
the objective function is no n-se parable, it still could be M-convex as long as 
the condition in Definition EH is satisfied. It should be also noted that SDA 
can apply to non-separable M-convex minimization problems while DA can 
not. 

8 Some algorithms in Q3-Q9) are based on the properties other than 
the optimality criterion, e.g., algorithm in ED is devised by utilizing the 
proximity theorem of M-convexity. The works in m also show how to 
obtain the starting point r(°) by constructing a convex extension of the M- 
convex objective function. 


Then, we have 


rN+ 

jex jey 


(n n (n 

jex jey 


(n u (n 

jex jey 


< 


n 


Hi 


jexuy 


(n ^ 

jexny 


( 20 ) 


which proves Lemma lA. 1 1 ■ 

Corollary A. 2: For Hj C V and y C /C such that X D 
y + 0 and X u y + JC, 


n«; + 

jeic\x 


n«; 

jeic\y 


s I n «? 

je{K\x)u(K\y) 


n «; 

je(K\x)n(K\y) 


Proof: The proof can be show by substituting X by K\X, 
y by K \ y and Hj by H'j in Lemma IA.1I ■ 
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